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ABSTRACT 


GOES-R series spacecraft feature a number of flexible appendages with modal frequencies 
below 3.0 Hz which, if excited by spacecraft disturbances, can be sources of undesirable jitter 
perturbing spacecraft pointing. In order to meet GOES-R pointing stability requirements, the 
spacecraft flight software implements an Active Vibration Damping (AVD) rate control law 
which acts in parallel with the nadir point attitude control law. The AVD controller commands 
spacecraft reaction wheel actuators based upon Inertial Measurement Unit (IMU) inputs to 
provide additional damping for spacecraft structural modes below 3.0 Hz which vary with 
solar wing angle. A GOES-R spacecraft dynamics and attitude control system identified 
model is constructed from pseudo-random reaction wheel torque commands and IMU angular 
rate response measurements occurring over a single orbit during spacecraft post-deployment 
activities. The identified Fourier model is computed on the ground, uplinked to the spacecraft 
flight computer, and the AVD controller filter coefficients are periodically computed on-board 
from the Fourier model. Consequently, the AVD controller formulation is based not upon 
pre-launch simulation model estimates but upon on-orbit nadir point attitude control and time- 
varying spacecraft dynamics. GOES-R high-fidelity time domain simulation results herein 
demonstrate the accuracy of the AVD identified Fourier model relative to the pre-launch 
spacecraft dynamics and control truth model. The AVD controller on-board the GOES-16 
spacecraft achieves more than a ten-fold increase in structural mode damping of the 
fundamental solar wing mode while maintaining controller stability margins and ensuring that 
the nadir point attitude control bandwidth does not fall below 0.02 Hz. On-orbit GOES-16 
spacecraft appendage modal frequencies and damping ratios are quantified based upon the 
AVD system identification, and the increase in modal damping provided by the AVD 
controller for each structural mode is presented. The GOES-16 spacecraft AVD controller 
frequency domain stability margins and nadir point attitude control bandwidth are presented 
along with on-orbit time domain disturbance response performance. 
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1 INTRODUCTION 


The Geostationary Operational Environmental Satellite-R (GOES-R) spacecraft was launched on 
November 19, 2016, and the spacecraft is now on station in geosynchronous orbit and renamed 
GOES-16. The Lockheed Martin GOES-R series spacecraft are the next generation of advanced 
geosynchronous weather satellites for the United States operated by the National Oceanic and 
Atmospheric Administration (NOAA). When fully deployed, the GOES-R series spacecraft feature 
a number of flexible appendages with modal frequencies below 3.0 Hz which, if excited by spacecraft 
disturbances, can be sources of undesirable jitter perturbing spacecraft pointing. In order to meet the 
GOES-R Guidance Navigation and Control (GN&C) subsystem pointing stability requirement of 221 
urad (or 46 arcseconds) peak-to-peak per axis over 60 seconds, the GOES-R flight software 
implements an Active Vibration Damping (AVD) rate control law which acts in parallel with the 
nominal mission Single Input Single Output (SISO) nadir point attitude control law. The AVD 
controller commands spacecraft reaction wheel actuators based upon Inertial Measurement Unit 
(IMU) inputs to provide additional damping for spacecraft structural modes below 3.0 Hz which vary 
with solar wing angle. 


A GOES-16 spacecraft dynamics and attitude control system identified model is constructed from 
pseudo-random reaction wheel torque commands and IMU angular rate response measurements 
occurring over a single orbit during spacecraft post-deployment activities. Ground processing 
algorithms perform a least-squares Fourier fit of these data as a function of solar wing angle to identify 
the exact character of the combined spacecraft on-orbit structural dynamics and nadir point attitude 
control law. The Fourier model coefficients are uplinked and stored on-board the spacecraft flight 
computer and used over the entire mission duration; at every two-degrees of solar wing rotation, the 
AVD controller filter coefficients are updated on-board from the identified Fourier model. 
Consequently, the AVD controller formulation is based not upon pre-launch simulation model 
estimates but upon on-orbit nadir point attitude control and time-varying spacecraft dynamics. 


This document includes four main sections, numbered as Section 2 through Section 5. Section 2 
presents an overview of the GOES-R series spacecraft deployed configuration which identifies the 
flexible appendage structures. Section 3 describes the three component parts of the AVD algorithm: 
on-orbit AVD excitation, ground-based AVD system identification, and on-orbit AVD compensation. 


GOES-R series spacecraft high-fidelity time domain simulation results in Section 4 demonstrate the 
accuracy of the AVD identified Fourier model relative to the pre-launch spacecraft dynamics and 
control truth model. High-fidelity time domain simulation results also demonstrate the robust 
performance of the AVD controller in the presence of worst-case appendage modal parameter 
variations anticipated over the GOES-R spacecraft 15-year mission duration. 


Section 5 presents extensive frequency and time domain on-orbit performance of the AVD controller 
for the GOES-16 spacecraft. The AVD identified Fourier model for the GOES-16 spacecraft is used 
to predict on-orbit disturbance rejection transfer functions, on-orbit modal parameters, on-orbit AVD 
controller stability margins, on-orbit nadir point attitude control bandwidth, and on-orbit impulse 
disturbance responses in the time domain. Finally, GOES-16 attitude error flight telemetry during 
multiple thruster firings demonstrates the on-orbit appendage vibration attenuation performance of 
the AVD controller. 
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2 GOES-R SERIES SPACECRAFT CONFIGURATION 


The fully deployed GOES-R series spacecraft configuration is shown in Figure | and Figure 2 along 
with the spacecraft body axis coordinate frame. Figure 2 documents the positive solar wing rotation 
angle convention, the solar wing zero-angle position, and the sun-facing side of the solar array. 
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Figure 1. GOES-R Series Spacecraft in the Deployed Configuration 
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Figure 2. Deployed GOES-R Series Spacecraft and Solar Wing Angle Convention 
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Each spacecraft carries a number of nadir pointing and solar pointing instruments as its primary 
payload. Each GOES-R series spacecraft features flexible appendages having bending modes below 
3.0 Hz when fully deployed. The AVD system identification process identifies five solar wing natural 
modes, two magnetometer boom modes, and two antenna wing modes. The AVD identified Fourier 
model is used to quantify modal parameters for these nine flexible appendage modes, and the AVD 
controller provides additional damping to attenuate vibrations for the majority of these appendage 
natural modes. 


3. ACTIVE VIBRATION DAMPING ALGORITHM 


The AVD controller designed for the GOES-R series spacecraft is a system-identification based high- 
bandwidth angular rate controller which acts in parallel with the attitude controller. AVD improves 
instrument pointing performance by applying phase stabilizing gain at the natural frequencies of the 
identified spacecraft appendage modes which modifies the spacecraft system dynamics, attenuates 
appendage vibration, and increases modal damping of the flexible appendages. AVD gain stabilizes 
all other spacecraft appendage modes with adequate margin to satisfy control system stability 
requirements. AVD also provides additional low-frequency disturbance attenuation control for all 
spacecraft axes which improves the performance of the nadir point attitude control law. AVD for the 
GOES-R series spacecraft is implemented as a single-input single- 
On-Orbit AVD output (SISO) rate controller used for each vehicle axis separately. 
i LL The AVD algorithm consists of three components which will be 
discussed separately in Section 3.1 through Section 3.3: a spacecraft 
Flight Telemetry flight software algorithm performing on-orbit AVD excitation, a 
RWA Torque at 20 Hz ‘ ieee ae 
IMU Rate at 20 Hz ground-based AVD system identification process, and a spacecraft 
SADA Angle at 1 Hz flight software algorithm performing on-orbit AVD compensation, as 
shown in Figure 3. GOES-R spacecraft flight telemetry from the 
Ground-Based AVD excitation process is passed to the ground-based AVD system 
AVD System identification process where the AVD identified Fourier model is 
Identification computed as a function of solar wing angle. Additional ground 
processing is performed using the Fourier model to select the AVD 
controller gains. The AVD identified Fourier model coefficients and 
the AVD controller gains are subsequently uplinked to the spacecraft 
and stored in flight computer memory. The AVD compensation 
flight software algorithm periodically updates the AVD controller 
i filter coefficients, which are computed from the identified Fourier 
Sabet model coefficients, the solar wing angle, and the AVD controller 
gains. 


Uplink to Flight Computer 
AVD Controller Gains 
Fourier Model Coefficients 


Compensation 


Figure 3. AVD Algorithm 
Components 


The AVD compensation algorithm, used on-board the spacecraft flight computer to compute the AVD 
controller filter coefficients, is used during ground processing to characterize all aspects of AVD 
controller performance after the identified Fourier model coefficients are computed. AVD controller 
disturbance rejection, stability margins, nadir point attitude control bandwidth, and time domain 
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impulse responses are analyzed prior to uplinking the Fourier model coefficients and AVD controller 
gains to the spacecraft flight computer. Samples of these results for the GOES-16 spacecraft are 
documented in Section 5. 


The AVD controller for the GOES-R series spacecraft is illustrated in the Figure 4 block diagram. 
When the AVD controller is disabled, the nominal nadir point controller provides attitude and rate 
control to stabilize spacecraft pointing. Upon enabling the AVD controller, the AVD compensation 
torque augments the nadir point attitude control torque, and both controllers provide combined 
attitude and rate control to stabilize the spacecraft and attenuate appendage vibration. 


t) 
O) Attitude OC} T + 7 Spacecraft 
C) Control Law ie. Dynamics 


Figure 4. AVD Rate Controller and Closed-Loop Spacecraft Dynamics 


3.1. On-Orbit AVD Excitation Process 

GOES-R spacecraft dynamics and control are identified once during the spacecraft post-launch on- 
orbit testing period and used over the entire mission duration. Because the GOES-R nadir point 
attitude control law is active during the system identification test, the closed-loop spacecraft dynamics 
and nadir point attitude control are identified together as a single Fourier model. To characterize the 
dependence of the spacecraft dynamics upon the solar wing angle, the system dynamics are identified 
over one complete solar wing rotation. 


Figure 5. Spacecraft Configuration during AVD Excitation for System Identification 


The spacecraft dynamics are excited by pseudo-random torques computed on-board and commanded 
to Reaction Wheel Assemblies (RWAs) during the on-orbit identification testing, as illustrated in 
Figure 5. AVD excitation is initiated by ground command. The response of the spacecraft to these 
inputs is determined by the IMU, which measures spacecraft body axis angular rates. The torque 
inputs and angular rate outputs telemetered at a 20 Hz sample rate are post-processed on the ground 
to identify the spacecraft dynamics as a function of solar wing angle. Identifying the spacecraft on- 
orbit dynamics and control eliminates the dependence of the AVD controller performance on the 
accuracy of pre-flight system dynamic models because the dynamics of the spacecraft, sensors, 
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actuators, and control processing are measured directly. The AVD identified model is represented 
by the shaded portion of Figure 5, which describes the torque input, T;,, to angular rate output, Woy, 
system model for each spacecraft body axis. 


AVD excitation is performed for each spacecraft axis separately every 10° of solar wing rotation, and 
each excitation interval lasts five minutes. The torque excitation zero-to-peak magnitude used for 
GOES-16 was 0.9 Nm per spacecraft body axis, and a representative angular rate response for the 
GOES- 16 roll axis at a 90.5° solar wing angle is shown in Figure 6. AVD excitation telemetry used 
for system identification consists of 36 torque input and 36 angular rate output time records per 
spacecraft body axis. 


Roll Dynamic Model Input Time History at SA Angle (deg) = 90.5 ° 10° Roll Dynamic Model Output Time History at SA Angle (deg) = 90.5 
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Output Angular Rate (rad/sec) 
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Figure 6. AVD Excitation Torque Input and Angular Rate Response Output for GOES-16 


3.2. Ground-Based AVD System Identification Process 


The measured spacecraft body rates, corresponding excitation torques, and solar wing angle are 
telemetered to the ground during post-launch testing of the spacecraft. This data is processed to 
estimate SISO transfer functions from RWA torque command to IMU angular rate response for each 
spacecraft body axis. The input-output data is collected for one full orbit period in order to capture 
the variation of the spacecraft dynamics over a 360° range of solar wing angles. 


Spacecraft Dynamics and Attitude Control Identified by AVD Excitation 


Figure 7. The AVD Identified Fourier Model 
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The first objective of the system identification process is to use the telemetered spacecraft data to 
generate the AVD identified Fourier model illustrated in Figure 7. The identified model consists of 
three independent and uncoupled SISO system models, one for each spacecraft body axis, from 
spacecraft torque command input, T;,, to angular rate response output, Woyr. The shaded block 
diagram encompassing the closed-loop nadir point attitude control law and spacecraft system 
dynamics is used throughout this document to always indicate the AVD identified Fourier model 
generated by the system identification process described in this section. 


The discrete-time AVD identified model transfer function, G(z), for each spacecraft body axis is 
assumed to be of the form shown in Equation 1. It is a p“” order RWA excitation torque, T(z), to 
IMU angular rate response, w(z), frequency domain transfer function. This is equivalent to the p“” 
order ARX (auto-regressive with exogenous input) discrete time domain model in Equation 2. 


biz 7t4+boz7 2 4+:+dyzP 
Giz) - 22-4 ee eee (1) 


T(z) 1=a,27-4=a32-7="-=a,2P 
w(k) = ayw(k — 1) + agw(k — 2) +--+ dpw(k — p) + by T(k — 1) + boT(k — 2) +--+ + byT (kK — p) (2) 


The standard approach for identifying the parameters of an ARX model is extended in order to capture 
the dependence of spacecraft dynamics upon the solar wing angle, 8. Equation 3 extends the ARX 
formulation of Equation 2 for the GOES-R series spacecraft to include periodicity of the identified 
model coefficients with respect to the solar wing angle. Each identified model coefficient in Equation 
2 is represented by a ninth-order Fourier series in Equation 3. Sine and cosine terms follow the 
compact convention sN6@; = sin(N@;) and cN@, = cos(N@,). The k subscript of the solar wing 
angle, 0,, indicates the angular rate response at discrete time, k, also depends upon the solar wing 
angle at time k. 


wk) = [Ayo + A11SOx% + A12COx + A13S20% + Ay4C20% + AysS30% + A16C30,% + A17S40x% + Aigc40,]a(k — 1) 
+ [Ago + Any SO, + Az2CO,% + A320, + Az4C20,% + Azs530% + Ar6C30,% + Az7S40,% + Azgc46;]w(k — 2) 
ae is 
+ [Apo + Ap1SOx + Ap2cOx + Ap3S2Ox + Apac2O% + Apss3O% + Ap6C3Ox + Ap7S40x + Apsc4Oxlo(k -—p) (3) 
+ [By + By SO, + By2C0;, + By3520; + By4c20;, + By5s30, + By6C30, + By7S40; + Bygc40;,|T(k — 1) 
+ [Bz + Bz1SO_ + Bz2CO, + Bz3S20;, + Bz4C20; + Bz5530, + BygC30,% + Bz7S40;, + Bzgc40;,|T(k — 2) 
fp 33 
+ [Boo + By1SOx + By2CO, + By3820, + By4C2O, + BysS30, + BygC30, + By7S40; + Bygc4O;,|T(k — p) 


Equation 3 is the formulation of the SISO discrete-time AVD identified Fourier model for each 
spacecraft body axis. The GOES-R series spacecraft use a 140" order model, p = 140, which yields 
1260 Aj; and 1260 B;; coefficients per spacecraft body axis. The total number of Fourier coefficients 
is 7560 encompassing all three spacecraft axes. As mentioned earlier, these coefficients are uplinked 
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to the spacecraft flight computer and stored in flight computer memory. The Fourier coefficients are 
computed by a standard least-squares solution approach performed during ground processing. 
Equation 3 is re-organized into the form shown in Equation 4 in order to separate the Fourier 
coefficients from the time dependent terms. Present values of the angular rate response, w(i), are 
ordered on the left side of the equation in column matrix Q, and all past time history values from 
flight telemetry are multiplied by sine and cosine multiples of the solar wing angle, 6;, and grouped 
into the huge matrix, O. The least-squares solution of the AVD identified Fourier model coefficients, 
®, is performed according to Equation 5. 


Q. =|] = O@ = |w(i- 1) s6,w(i-1) c6,w(i-—1) 520,00 — 1) c20,w(i—1) ... s40,TG— p) c40,T(i— p)]| (4) 


® = [0’0]-‘0'70 (5) 


After computing all Fourier coefficients of the identified model in Equation 3, the practical use of the 
equation becomes clear. Given a particular GOES-R solar wing angle, 0, the Fourier series 
coefficients in Equation 3 reduce to single a; and b; coefficients of Equation | and Equation 2, which 
can then be used for frequency domain or time domain analysis, respectively, of the identified system 
dynamics and control. 


As mentioned earlier, the AVD controller filter coefficients are a function of the identified Fourier 
model coefficients, the solar wing angle, and the AVD controller gains. The development of the 
AVD control law formulation is proprietary technical information owned by Lockheed Martin and 
patented by two of the authors [1]. The algorithm to compute the AVD controller filter equations is 
implemented in the GOES-R series spacecraft flight computer. The algorithm computes the 
coefficients of the controller transfer function given by Equation 6, and the AVD control torque is 
computed using the discrete-time representation of Equation 7. 


The SISO AVD controller is also a set of three 140" order discrete transfer functions, one for each 
spacecraft body axis. The AVD controller coefficients g; and h; are a function of GOES-R solar 
wing angle, as they are computed directly from the AVD identified model coefficients a; and b; in 
Equation 2, which are, in turn, computed from Equation 3 given the solar wing angle of interest. The 
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compensation provided by the AVD controller uniquely applies phase stabilizing gain at the natural 
frequencies of the identified spacecraft appendage modes to attenuate appendage vibration; all other 
spacecraft structural modes are gain stabilized. 


_ 7) _ G12 * +9227? +-+Gp2? 
CO aaa ee 6) 
T(k) = hy T(k ~ 1) + hgT(k = 2) + + pT (k ~ p) + gra(k ~ 1) + gawk — 2) ++ gpalk ~p) (7) 


3.3. On-Orbit AVD Compensation Process 


After completing system identification and verifying the performance and robustness of the AVD 
identification-based angular rate controller, the identified Fourier model is stored on-board the 
GOES-R series spacecraft flight computer. The identified model coefficients are used on-board to 
compute the AVD controller filter coefficients, which are updated every two degrees of solar wing 
rotation (every 8 minutes). The AVD controller receives IMU input at 20 Hz and computes 20 Hz 
control torque commands, which are summed with other feed-forward torque commands generated 
by the GOES-R attitude control system. AWD compensation is always active during nadir point 
operation of the GOES-R series spacecraft. 


When the 140" order AVD controller was first activated on-board the GOES-16 spacecraft on 
February 1, 2017, a perceptible 1.1% increase in computational load was observed for the RAD750 
flight computer processor. The computational load increase occurred because the AVD compensation 
torque per spacecraft body axis (Equation 7) is computed at 20 Hz from 140 past values of AVD 
torque and 140 past values of IMU angular rate measurements. 


4 AVD PRE-FLIGHT HIGH-FIDELITY SIMULATION FOR GOES-R 


The Lockheed Martin guidance and control team for the GOES-R program developed a high-fidelity 
spacecraft pointing control and jitter time domain simulation used to verify GOES-R spacecraft 
requirements during the design and development phase of the program. The high-fidelity simulation 
includes all aspects of the pointing control subsystem, and the simulation includes a detailed structural 
dynamic model of the GOES-R series spacecraft used to characterize the attitude control and 
structural dynamic interaction predicted for flight. 


4.1 AVD Identified Model Transfer Functions versus Spacecraft Dynamics Open-Loop Truth Model 


The on-orbit AVD identified Fourier model for GOES-R series spacecraft is one of the output 
products derived from the AVD system identification process. The identified model represents the 
true characteristics of the on-orbit spacecraft dynamics and control, but how accurately does the 
identified model predict the modal character of the spacecraft appendages? The AVD algorithm was 
implemented in the high-fidelity simulation to perform AVD excitation as described in Section 3.1, 
and the AVD ground processing system identification tools (Section 3.2) were used to post-process 
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the simulated AVD excitation data to generate an AVD identified Fourier model. The identified 
model is compared with the open-loop system dynamic model from the high-fidelity simulation to 
demonstrate the ability of the AVD identified Fourier model to correctly characterize system 
appendage mode dynamics of the truth model. Examples of disturbance response transfer functions 
from body frame torque input to angular rate response output are shown in Figure 8 for the roll and 
yaw axes at a 0° solar wing angle. As expected, the AVD identified model transfer functions nearly 
exhibit the characteristic -20 dB per decade low-frequency roll-off associated with the attitude 
controller angular rate response, but the key characteristic illustrated in these figures is the ability of 
the AVD identified Fourier model to accurately predict the appendage modal peaks. 


‘ Roll Axis Original vs. Identified TF at SA Angle (deg) = 0.0 . Yaw Axis Original vs. Identified TF at SA Angle (deg) = 0.0 
BO i ng a a 
\_— Open Loop Original Model | E ——— Open Loop Original Mode! 


f ——— Open Loop identified Mode! } f ——— Open Loop Identified Mode! 


Angular Rate vs. Torque Magnitude (rad/sec)/(N*m) 
Angular Rate vs. Torque Magnitude (rad/sec)/(N*m) 


7| Discrete Transfer Function ID Model Order = 200 L 7| Discrete Transfer Function ID Model Order = 200 LL i 


10 & : F ar 2 ' 
10 10° 10 10 10 10 10 10 
Frequency (Hz) Frequency (Hz) 


Figure 8. Open-Loop Spacecraft Dynamic Truth Model vs. AVD ID Model at 0° SADA Angle 


0 


4.2 AVD Identified Model Modal Parameters versus Spacecraft Dynamics Open-Loop Truth Model 


Consistent with the disturbance response transfer functions presented in Section 4.1, the modal 
parameters are tabulated for nine appendage modes of the GOES-R series spacecraft. The modal 
parameters of the simulated AVD identified Fourier model are compared with the high-fidelity 
simulation open-loop system dynamics truth model. 


Table 1: Modal Parameter Comparison Open-Loop Dynamics Truth vs. AVD ID Model 


Dynamics Truth Model || AVD Identified Model || Percent Difference 


Appendage Flexible Mode Description .. ||Frequency| Damping | Frequency | Damping ee Damping 
(Hz) | Ratio(%) | (Hz) | Ratio(%) | OY 


Solar Wing Out-of-Plane Bending | | o | 
Magnetometer Boom X-Rotation Bending | 330. | 


Magnetometer Boom Y-Rotation Bending 
Solar Wing Out-of-Plane Bending II 


All 

| SolarWingin-PlaneBending | 0 
| Solarwingtwist 
| Solar Wing Out-of- Plane Bending | 0 
| Antenna Wingin-PlaneBending |_| 
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Table 1 lists the modal frequencies and damping ratios for nine appendage modes of the high-fidelity 
simulation structural dynamics truth model and the AVD identified model. Modal parameter 
percentage differences are also tabulated. The AVD identified model determines the appendage 
modal frequencies within a 1.2% error; however, the identified model predicts modal damping within 
a 22% error. Because the AVD identified model is used in the AVD rate controller formulation to 
apply phase stabilizing gain at the spacecraft appendage modal frequencies, AVD controller 
effectiveness depends more upon accurate knowledge of appendage modal frequencies than damping 
ratios. 


4.3 AVD Controller Robustness Analysis 


Robustness of the AVD controller to appendage modal frequency variation was quantified using the 
GOES-R high-fidelity pointing control and jitter time domain simulation. Over the GOES-R series 
spacecraft 15-year mission duration, no more than a 5% variation in appendage modal frequencies is 
anticipated due primarily to spacecraft mass and inertia changes as thruster propellant is depleted. 
Because the AVD identified Fourier model is derived once during post-launch testing and used over 
the entire mission duration, AVD controller effectiveness may decline as spacecraft appendage modal 
frequencies change over time. 


Quantifying AVD controller robustness required generating 160 high-fidelity simulation structural 
dynamic models with uniformly varying appendage modal frequencies within +5% of the pre-launch 
nominal values. All simulations used the same AVD identified Fourier model to generate AVD 
controller filter coefficients based upon the nominal structural dynamic model. Time domain 
simulation scenarios evaluated GOES-R spacecraft solar instrument pointing errors, solar instrument 
pointing stability, nadir instrument peak attitude and angular rate errors, and nadir instrument pointing 
stability. Simulation study results demonstrated robust AVD controller performance in all cases, all 
simulation results remained compliant with all GOES-R pointing requirements, and spacecraft 
pointing performance was only slightly affected by appendage modal frequency variation with the 
nominal AVD controller. Even though the AVD excitation test is planned only once during post- 
launch calibrations, AVD excitation calibration can be repeated at any time during the spacecraft 
mission to generate a new AVD identified Fourier model if spacecraft disturbance time responses 
begin to show unacceptable appendage vibration decay times affecting instrument pointing 
performance. 


5 AVD ON-ORBIT FLIGHT PERFORMANCE FOR GOES-16 


GOES-16 on-orbit AVD excitation was performed on 5-7 January 2017. Ground-based AVD system 
identification analysis commenced thereafter, and the AVD identified Fourier model coefficients and 
AVD controller gains were uplinked to the GOES-16 spacecraft on 1 February 2017. All frequency 
and time domain results presented in this section are based entirely upon the AVD identified Fourier 
model generated from GOES-16 AVD excitation flight telemetry. Section 5.1 through Section 5.5 
present predicted on-orbit performance using the GOES-16 AVD identified model, whereas on-orbit 
AVD controller time domain performance from GOES-16 flight telemetry is presented in Section 5.6. 
These results demonstrate the spectacular on-orbit AVD controller performance for the GOES-16 
spacecraft. 
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5.1 GOES-16 System-Identified Model Transfer Functions with AVD Disabled vs. AVD Enabled 


Disturbance torque input, T;,y, to angular rate response output, Woyr, frequency domain transfer 
functions are generated for the GOES-16 AVD identified model based upon Figure 7; these are the 
GOES- 16 on-orbit disturbance response transfer functions in nadir point control mode with the AVD 
controller disabled. For a particular solar wing angle, the GOES-16 AVD controller filter coefficients 
are computed by Equation (7) in Section 3.2, and the AVD controller loop is closed around the AVD 
identified model as shown in Figure 9. Disturbance torque input, T¢,, to angular rate response output, 
Wour, frequency domain transfer functions, computed using Figure 9, are the on-orbit disturbance 
responses of the GOES-16 spacecraft in nadir point control mode with the AVD controller enabled. 
Figure 10 and Figure 11 plot the GOES-16 disturbance rejection transfer functions which illustrate 
the increase in appendage modal damping with the AVD controller enabled as compared with the 
AVD controller disabled. 


Spacecraft Dynamics and Attitude Control Identified by AVD Excitation 


Figure 9. Closed-Loop AVD Controller with AVD Identified Model 
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Figure 10. Roll Axis Disturbance Rejection with AVD Disabled vs. Enabled 


Many of the GOES-16 disturbance torque input to angular rate output transfer functions generated 
using the AVD identified model do not show the ideal -20 dB per decade low-frequency roll-off down 
to 0.001 Hz, as illustrated in Figure 10 and Figure 11. Because 5-minute AVD excitation intervals 
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are used for the GOES-16 AVD calibration, the minimum frequency content present in the AVD 
excitation telemetry data is 0.0033 Hz. The low-frequency disturbance response characteristic is not 
a pointing performance concern because the AVD controller acts in parallel with the nadir point 
attitude control law, as shown in Figure 9, and the nadir point attitude controller attenuates very low 
frequency disturbances. 
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Figure 11. Pitch and Yaw Axis Disturbance Rejection with AVD Disabled vs. Enabled 


5.2. GOES-16 System-Identified Model Modal Parameters with AVD Disabled vs. AVD Enabled 


GOES- 16 appendage modal parameters associated with the transfer functions shown in Figure 10 and 
Figure 11 are listed in Table 2. Modal parameters with the AVD controller disabled are the modal 
parameters of the AVD identified model. The modal parameters in Table 2 are representative of the 
AVD controller performance for the GOES-16 spacecraft for all solar wing angles. 


Table 2: GOES-16 Modal Parameters AVD Disabled vs. AVD Enabled 


eae Damping | Damping AVD 
Appendage Flexible Mode Description me gerana! Ratio (%) | Ratio(%) | Damping 
Axis (Hz) oe ‘i 
Identified | with AVD Factor 


[Solar wing Out-of Plane Bending) [0 | Roll | 027 | 0253 | 2703 | 1082 
[ Magnetometer 80om x-Rotation Bending | 330 | Roll | o6s1 | om | 0.628 | 0.93 
[magnetometer Boom Y-Rotation Bending | 250 | itch | 0.733 | 0.820 | 0.926 | 1.01 
Solar Wing Out-of Plane Bending | 0 
antenna Wing Out-of Plane Bending | all | pitch | 0999 | 7191 | 782 | 1.00 
[solar wingin-Plane Bending [0 | vaw | 1.313 | 0978 | 1971 | 202 
FsolarwingTwist __-—~—~—S~w CT itch | 1.320 | 0279 | 0366 [131 
Solar Wing Out-of-Plane Bending [0 | roll | 1923 | 0321 | 0305 | 0.95 
fo | Yew | 1.996 | 040 | 1149 [1.37 


Antenna Wing In-Plane Bending 


Note: AVD Damping Factor = (Modal Damping AVD Enabled)/(Modal Damping AVD Disabled) 


The AVD controller does not increase modal damping for all appendage modes, but five of the nine 
appendage structural modes identified by AVD excitation realize an increase in modal damping. In 
particular, the AVD controller increases the modal damping for the GOES-16 fundamental solar wing 
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mode (the first out-of-plane bending mode of the solar wing) by more than a factor of 10. This is the 
benefit provided by AVD; the AVD controller is a flight software algorithm effectively adding 
damping to GOES-16 appendage hardware. 


On-orbit GOES-16 appendage modal frequencies identified by AVD excitation are listed in Table 2. 
These appendage modal frequencies are higher than the pre-launch model-predicted appendage 
frequencies for GOES-R series spacecraft listed in Table 1. The lower model-predicted frequencies 
are due to conservatism built-in to structural dynamic finite element models developed during the 
GOES-R spacecraft design phase. This comparison illustrates a key benefit of the AVD controller 
formulation based upon on-orbit system identification dynamics and control models rather than less 
accurate pre-launch simulation models. 


5.3. GOES-16 Frequency Domain AVD Controller Stability Margins 


GOES-16 AVD controller stability margins are determined from the open-loop transfer function of 
the AVD controller in series with the AVD identified Fourier model, as illustrated by the block 
diagram in Figure 12. 


Spacecraft Dynamics and Attitude Control Identified by AVD Excitation 


Figure 12. AVD Controller in Series with AVD Identified Model 


Frequency domain transfer function Nichols plots from W ;y to Woyr plotted as a magnitude response 
in decibels versus phase angle in degrees are shown in Figure 13. These figures are representative 
Nichols plots generated to characterize the GOES-16 AVD controller at all solar wing angles at 1° 
increments for all three spacecraft body axes. Figure 13 shows the frequency domain characteristics 
of the AVD controller as it applies phase stabilizing gain at the GOES-16 appendage natural 
frequencies to quickly attenuate appendage vibration measured by the spacecraft IMU. Phase 
margins in degrees are measured along the zero decibel line between the transfer function modal 
peaks and the vertical lines at 360° multiples of the -180° line. Phase margins are listed on each 
figure. Gain margins are the gain in decibels from the zero decibel line to the transfer function curves 
as they cross vertical lines at 360° multiples of the -180° line. AVD controller minimum gain and 
phase margins at each solar wing angle at 1° increments are evaluated and plotted versus the solar 
wing angle in Figure 14 for all three GOES-16 spacecraft body axes. The smallest gain margin is 
11.5 decibels in the roll axis at a 41° solar wing angle, and the smallest phase margin is 78° also in 
the roll axis at a 219° solar wing angle. The large stability margins for the GOES-16 AVD controller 
are another indication of the robustness of the vibration controller design. 
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Figure 13. GOES-16 AVD Controller Stability via Roll and Yaw Axis Nichols Plots 


Minimum Gain Margins vs. SADA Angle Minimum Phase Margins vs. SADA Angle 
20 1 T T T T 180 T 1 T T — i 
Roll Axis / Roll Axis 
rs Pitch Axis Pitch Axis 
—— Yaw Axis 1607 Yaw Axis 
Gain Margin Limit | Phase Margin ul 
10 es 
140+ | 
5 4 
11.5 dB Minimum Gain Margin @ 120+ | Fy 
a Roll Axis at 41° Solar Wing Angle 3 F 
@ aL | 2 \ oY 
$0 g et g 
< Q 100+ \ / 
>, £ % pa 
\ / 
= 2 ot SY 
c | Ss 80 
BS -10/ | ® 
co) -— 8 
/ a EOS yee OS ~< 607 78° Minimum Phase Margin in Roll Axis at 219° Solar Wing Angle 
15ST 6s . Al . a 
5 L/ gis a i NY 7 \ =< = for 0.303 Hz Solar Wing Out-of-Plane Fundamental Mode 
x xX \/ il Sf 
hae \ VW : L~m -—7 40+ 
-20P~ \ ~ fo | 
Wieser \ sale 
waa / Xp | 
-25} OL System: AVD Control in Serids with ID Closed-Loop Attitude Control | 207 on s cee AVD pears in Series with ID Closed-Loop Attitude Control 
AVD Gamma Gain = [ 0.9811, 0.9839, 0. 9731 ] AVD ain = [ 0.9811, 0.9839, 0.9731 ] 
AVD Controller Output Delay (sec) = 0.05 AVD Controller Output “Belay (sec) = 0.050 
30 Discrete Transfer Function ID Model Order = [ 140, 140, 140 ] i 0 Discrete Transfer Function ID Model Order = [ 140, 140, 140 ] i 
it) 50 100 150 200 250 300 350 0 50 100 150 200 250 300 350 
SADA Angle (Degrees) SADA Angle (Degrees) 


Figure 14. GOES-16 AVD Controller Minimum Gain and Phase Margins vs. SADA Angle 


5.4 GOES-16 Nadir Point Attitude Control System Bandwidth with AVD Disabled vs. AVD Enabled 


The AVD angular rate controller flattens the GOES-16 disturbance response curve as it provides 
additional disturbance rejection at the attitude control system bandwidth between 0.03-0.04 Hz, as is 
evident in Figure 10 and Figure 11. This effect also influences the nadir point attitude control system 
bandwidth with the AVD controller enabled. 


The bandwidth of the GOES-16 nadir point attitude controller is computed from the commanded 
attitude, 0-, to attitude response, 0, transfer function as illustrated in the Figure 15 block diagram and 
written as Equation 8. This diagram divides the attitude controller into its constituent parts: the 
attitude control path, G,, and the rate control path, Gp. The spacecraft dynamics is represented by the 
D transfer function block. The familiar AVD identified Fourier model block diagram using the same 
convention is shown in Figure 16. The AVD identified model transfer function from disturbance 
torque input, T;,, to angular rate response output, Woyr, is written as Equation 9. Equation 10 shows 
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that the transfer function from commanded attitude, 0¢, to attitude response, 8, can be constructed by 
augmenting the AVD identified model with the negative derivative of the attitude controller attitude 
path transfer function, (— G,/s). Therefore, the AVD identified model can be used to estimate the 
nadir point attitude control system bandwidth by augmenting it with an analytical expression for 
(— G,/s), as shown in block diagram form in Figure 17. 


Spacecraft Dynamics and Attitude Control Identified by AVD Excitation 


Figure 16. GOES-16 AVD Identified Fourier Model with Attitude Control Law Divided into 
Constituent Parts 
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Spacecraft Dynamics and Attitude Control Identified by AVD Excitation 


Figure 17. GOES-16 AVD Identified Model Augmented to Estimate Nadir Point Attitude 
Control System Bandwidth with AVD Compensation Disabled 


Figure 17 is used to compute the commanded attitude, 9c, to attitude response, 0, transfer function 
for the GOES-16 nadir point attitude control system with the AVD controller disabled, and Figure 18 
is used to compute the same transfer function with the AVD controller enabled. Representative 
examples these transfer functions for GOES-16 are shown in Figure 19. The roll axis nadir point 
attitude control system bandwidth is 0.0444 Hz at the -3dB response amplitude for a 49° solar wing 
angle with the AVD controller disabled; with the AVD controller enabled, the bandwidth drops to 
0.02 Hz. A similar effect is observed for the yaw axis at a 22° solar wing angle in Figure 19. Because 
the GOES-16 primary mission operations concept is to stare at Earth while in nadir point attitude 
control mode, a drop in attitude control system bandwidth does not affect mission performance. The 
GOES-16 spacecraft is not required to quickly slew from one target to the next, which would require 
a high-bandwidth attitude control system. 


Spacecraft Dynamics and Attitude Control Identified by AVD Excitation 


Figure 18. GOES-16 AVD Identified Model Augmented to Estimate Nadir Point Attitude 
Control System Bandwidth with AVD Compensation Enabled 
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Figure 19. GOES-16 Nadir Point Attitude Control System Bandwidth with AVD 
Compensation Disabled and Enabled for Roll and Yaw Axes 


Enabling the AVD controller always causes a drop in nadir point attitude control system bandwidth 
compared with the attitude control system bandwidth when AVD is disabled. Augmenting the AVD 
identified model with the (— G,/s) transfer function effectively rotates the GOES-16 disturbance 
rejection curves (Figure 10 and Figure 11) clockwise and applies gain to shift the rotated disturbance 
rejection curves up to unity magnitude. Therefore, as the AVD controller increases disturbance 
rejection at the attitude control system bandwidth (0.03-0.04 Hz), the AVD controller likewise causes 
a decrease in the nadir point attitude control system bandwidth, as Figure 19 shows. Also, because 
the AVD identified model angular rate disturbance rejection curves do not follow the typical -20 
decibels per decade low-frequency roll-off, the attitude command to attitude response transfer 
functions in Figure 19 are not flat in the low frequency portion of the spectrum. 
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Figure 20. GOES-16 Nadir Point Attitude Control System Bandwidth with AVD 
Compensation Disabled and Enabled for All Solar Wing Angles 
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Using the method described in this section, the GOES-16 spacecraft nadir point attitude control 
system bandwidth over all solar wing angles at 1° increments with the AVD controller disabled and 
enabled is shown in Figure 20. The on-orbit nadir point attitude control system bandwidth with AVD 
disabled is approximately 0.045 Hz for all spacecraft body axes at all solar wing angles, whereas with 
AVD compensation enabled, the attitude control system bandwidth varies widely with solar wing 
angle. The AVD controller gains were selected to ensure that the nadir point attitude control system 
bandwidth is always greater than or equal to 0.02 Hz for all GOES-16 spacecraft body axes over all 
solar wing angles. 


5.5. GOES-16 Simulated Time Domain Impulse Responses with AVD Disabled vs. AVD Enabled 


Much of the analysis presented thus far was performed in the frequency domain using the frequency 
domain transfer functions for the GOES-16 AVD identified model (Equation 1) and the GOES-16 
AVD controller (Equation 6). This section presents time domain simulation results using the GOES- 
16 AVD identified model (Equation 2) and the GOES-16 AVD controller (Equation 7). Time domain 
simulation from disturbance torque input to angular rate response output using these GOES-16 
models follows the block diagram formats already presented when the AVD controller is disabled 
(Figure 7) or enabled (Figure 9). 


The time domain simulation using the GOES-16 flight models is forced by 0.05 second duration 
torque impulse inputs of 1.0 Nm amplitude for the roll and pitch axes and 0.5 Nm amplitude in yaw. 
The impulses occur every 24 minutes (every 6° of solar wing rotation), and the simulation time 
duration is one full day starting at a 0° solar wing angle. As occurs for GOES-16 in flight, the AVD 
identified Fourier model and the AVD controller filter coefficients are updated at every 2° of solar 
wing rotation. The simulated angular rate responses are shown in Figure 21 through Figure 23 for 
the roll, pitch, and yaw axes, respectively, with the AVD controller disabled and enabled. Each figure 
shows time responses for one full day along with expanded plots covering 124 minutes in duration. 
Overlaying the AVD enabled time responses on top of the AVD disabled responses illustrates the 
effectiveness of the AVD controller to quickly attenuate GOES-16 appendage vibrations over the 
simulated orbit. 
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Figure 21. Simulated GOES-16 Roll Angular Rate Responses due to Torque Impulses 
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Figure 22. Simulated GOES-16 Pitch Angular Rate Responses due to Torque Impulses 
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Figure 23. Simulated GOES-16 Yaw Angular Rate Responses due to Torque Impulses 


Figure 21 and Figure 23 show that roll and yaw axis angular rate responses due to torque impulse 
inputs can last tens of minutes in duration with the AVD controller disabled, but the AVD controller 
attenuates the appendage vibrations in approximately one minute or less. The AVD controller is 
effective in reducing pitch axis appendage vibrations (Figure 22), but the pitch axis vibration 
attenuation is not as dramatic when compared with the roll and yaw axes. Perhaps the most important 
results illustrated by Figure 21 through Figure 23 are that GOES-16 angular rate time responses with 
the AVD controller enabled are always bounded over all solar wing angles, and appendage vibration 
decay times with the AVD controller enabled are always less than with the AVD controller disabled. 


5.6 GOES-16 On-Orbit Time Domain Disturbance Response with AVD Disabled vs. AVD Enabled 

GOES- 16 on-orbit attitude error response is usually quiet and uneventful. However, every four days, 
arcjet thrusters fire for approximately one hour to perform north-south stationkeeping maneuvers, and 
low-thrust REA thrusters fire daily to unload reaction wheel momentum acquired during the previous 
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24 hours. These thruster firing events excite GOES-16 appendage vibration and provide a good 
opportunity to observe AVD controller performance. On 29 January 2017 and 17 February 2017, 
temporally similar thruster firing events occurred. A north-south stationkeeping maneuver finished 
30 minutes prior to the start of a momentum adjust maneuver. AVD compensation had not yet been 
enabled during the January 2017 event, but on the February 2017 occurrence, AVD compensation 
was enabled, which provided a means to compare GOES-16 appendage vibration responses as 
measured by the IMU. Figure 24 shows the GOES-16 roll axis attitude error response during these 
thruster firing events, and the pointing performance improvement with AVD compensation enabled 
is spectacular. Prior to activating AVD compensation, GOES-16 appendage vibration lasted tens of 
minutes after thruster firings, but with AVD compensation enabled, appendage vibration is quickly 
attenuated within one minute or less. The AVD algorithm is successfully operating on-board the 
GOES- 16 spacecraft. 
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Figure 24. GOES-16 On-Orbit Disturbance Response during Thruster Firings with AVD 
Compensation Disabled and Enabled 


6 CONCLUSION 


The Active Vibration Damping controller designed for the GOES-16 spacecraft is operational after 
successful implementation on-board the GOES-16 flight computer. The AVD flight software 
algorithm is an innovative control system formulation, generating reaction wheel torque 
compensation commands from IMU input, which achieves more than a ten-fold increase in damping 
for the GOES-16 solar wing fundamental appendage mode. The AVD controller designed for the 
GOES-16 spacecraft is a system-identification based high-bandwidth angular rate controller which 
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acts in parallel with the nadir point attitude controller. The AVD controller formulation is based not 
upon pre-launch simulation model estimates but upon on-orbit nadir point attitude control and time- 
varying spacecraft dynamics. AVD improves instrument pointing performance by applying robust 
phase stabilizing gain at the natural frequencies of the identified spacecraft appendage modes which 
modifies the spacecraft system dynamics, attenuates appendage vibration, and increases modal 
damping of the flexible appendages. AVD gain stabilizes all other spacecraft appendage modes with 
adequate margin to satisfy control system stability requirements. AVD provides additional low- 
frequency disturbance attenuation control for all spacecraft body axes which improves the 
performance of the nadir point attitude control law while ensuring that the attitude control bandwidth 
does not fall below 0.02 Hz. 
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